

%% specify files

matfiles = {'run1a_baseline_sr.mat','run1b_baseline_no_sr.mat',...
    'run5a_crra_sr.mat','run7a_Nstate_sr_fit_v3.mat','run6a_time_varying_s_sr.mat'};

for i = 1:5
    load(matfiles{i})
    stats_list{i} = stats;
end

%% moments
fprintf('\n\n')

fprintf('U: mean & %4.2f & %4.2f & %4.2f & %4.2f & %4.2f\n', cellfun(@(X)100*X.U_mean, stats_list));
fprintf('U: vol & %4.2f & %4.2f & %4.2f & %4.2f & %4.2f\n', cellfun(@(X)100*X.U_std, stats_list));
fprintf('Output gr vol& %4.2f & %4.2f & %4.2f & %4.2f & %4.2f\n', cellfun(@(X)100*X.Ygrowth_vol, stats_list));
fprintf('Cons gr ar1 & %4.2f & %4.2f & %4.2f & %4.2f & %4.2f\n', cellfun(@(X)X.Cgrowth_AR1, stats_list));
fprintf('Cons gr vol& %4.2f & %4.2f & %4.2f & %4.2f & %4.2f\n', cellfun(@(X)100*X.Cgrowth_vol, stats_list));
fprintf('1yr, mean & %4.2f & %4.2f & %4.2f & %4.2f & %4.2f\n', cellfun(@(X)X.nom_yld_mean(1), stats_list));
fprintf('1yr, std & %4.2f & %4.2f & %4.2f & %4.2f & %4.2f\n', cellfun(@(X)X.nom_yld_std(1), stats_list));
fprintf('wage vol/output vol & %4.2f & %4.2f & %4.2f & %4.2f & %4.2f\n', cellfun(@(X)X.stdwagebill2stdY, stats_list));
fprintf('corr(mu3) & %4.2f & %4.2f & %4.2f & %4.2f & %4.2f\n', cellfun(@(X)mean(X.cross_sec_cons.corr_mu3_1y_ny(1:5)), stats_list));
fprintf('vol(mu3) & %4.4f & %4.4f & %4.4f & %4.4f & %4.4f\n', cellfun(@(X)X.cross_sec_cons.std_mu3_1y, stats_list));
fprintf('Tightness, mean & %4.2f & %4.2f & %4.2f & %4.2f & %4.2f\n', cellfun(@(X)X.tightness_mean, stats_list));
fprintf('Tightness, std & %4.2f & %4.2f & %4.2f & %4.2f & %4.2f\n', cellfun(@(X)X.tightness_std, stats_list));
fprintf('corr(U,V) & %4.2f & %4.2f & %4.2f & %4.2f & %4.2f\n', cellfun(@(X)X.corr_U_V, stats_list));
fprintf('Eq prem, mean & %4.2f & %4.2f & %4.2f & %4.2f & %4.2f\n', cellfun(@(X)X.stock_ret_1mo.risk_prem_unlevered*3*12*100, stats_list));
fprintf('Eq prem, vol & %4.2f & %4.2f & %4.2f & %4.2f & %4.2f\n', cellfun(@(X)X.stock_ret_1mo.stock_exret_vol_unlevered*3*sqrt(12)*100, stats_list));
fprintf('DP, mean & %4.3f & %4.3f & %4.3f & %4.3f & %4.3f\n', cellfun(@(X)X.stock_ret_1mo.DPratio_mean*12, stats_list));
fprintf('DP, vol & %4.3f & %4.3f & %4.3f & %4.3f & %4.3f\n', cellfun(@(X)X.stock_ret_1mo.DPratio_std*12*3, stats_list)); % apply leverage factor of 3